1 



On 
On 



Dirac operator as a random matrix and the quenched limit of QCD with 
chemical potential 
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The behavior of quenched QCD at nonzero chemical potential p has been a long-standing puzzle. An explicit 
solution is found using the random matrix approach to chiral symmetry breaking. At nonzero p the quenched 
QCD is not a simple n — > limit of a theory with n quarks: a naive 'replica trick' fails. A limit that leads to the 
quenched QCD is that of a theory with 2n quarks: n quarks with original action and n quarks with conjugate 
action. 
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1. The puzzle 

The main result of the present study is a res- 
olution of a puzzle which existed since the mid- 
dle of 80's and has been one of the major ob- 
stacles towards our understanding of the behav- 
ior of QCD at finite baryon density. The contra- 
diction is simple and acute. The quark chemical 
potential p in QCD shifts the energy of a given 
quantum state by 3/i times this state's net baryon 
number. A transition at T = should occur 
when p — p c w m p /3: the vacuum state without 
baryons is no longer the lowest energy state. It is 
believed that the chiral symmetry is restored in 
the new phase. Testing this prediction by Monte 
Carlo simulation is difficult since the determi- 
nant of the fermion matrix becomes complex at 
nonzero p. Nevertheless one can use quenched 
approximation in this case. Such Monte Carlo 
studies jjj give consistently the critical value of 
p c « 771^/2 — a clearly unphysical result, since 
the pion does not carry baryon number. 

The present study is motivated by a desire 
to clarify the problem using random matrix ap- 
proach. On an exactly solvable model of chiral 
symmetry breaking in QCD we demonstrate ex- 
plicitly a general phenomenon related to the na- 
ture of the quenched limit. Namely, one can view 
the quenched theory as a somewhat abstract limit 
n — > of QCD with variable number n of quark 
flavors. Then if the limit is smooth the quenched 
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theory is a reasonable approximation to the the- 
ory at finite n. In fact, this limit is not always 
smooth and can be different from the quenched 
theory as, e.g., in a theory of spin glasses. We 
show that the quenched limit is not smooth in 
QCD at finite p. Quenched theory is therefore 
not a good approximation to full QCD at finite 
p. This resolves the old contradiction. More- 
over, we can describe a theory to which quenched 
QCD is a good approximation. It is a theory 
(suggested in ||) in which each quark has a con- 
jugate partner, essentially, with opposite baryon 
number assignment. 

2. (■ijjijj) an d P 

We study analytical properties of the quark 
condensate (V'V'} — the- order parameter of the 
chiral symmetry breaking — as a function of the 
bare quark mass m which we allow to go into the 
complex plane: m = z — x + iy. Integrating 
out the quark degrees of freedom we can express 
{ipip) through the average of the inverse Dirac op- 
erator D: (tpijj) = ( Tr (z - D)' 1 ) = G. In the 
infinite volume the value of the (V ; V-'} i n t ne cn i _ 
ral limit z — > can be related to the density of 
small eigenvalues p(X) by the Banks-Casher for- 
mula: ('tpip) — 7rp(0). 

To generalize this well-known relation to our 
case we must consider the density per unit area 
of the complex plane: p(x,y). The relation be- 
tween p and G is straightforward: G(x, y) = 
J dx'dy'p(x', y')(z — z')^ 1 . Viewing p as a density 
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of charges and Re G, — Im G as components of 
the electric field strength G one easily inverts this 
relation: 
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The main point here is that p vanishes where G 
depends analytically on z, the bare quark mass. 
For example, when p = all eigenvalues of D lie 
on the imaginary axis. This means that G has a 
cut as on Fig. |l|a. The discontinuity across the 
cut is the signature of the chiral symmetry break- 
ing and is given precisely by the Banks-Casher 
formula. 

3. Random matrices 

Now we would like to calculate G and p. Since 
we are interested in the density of small eigen- 
values of the Dirac operator one can expect that 
some simplification can be made. Indeed, it is 
by now well understood that these small eigen- 
values are related to the zero modes of the Dirac 
operator in the instanton field background. The 
behavior of these small eigenvalues is rather uni- 
versal and can be described by a simple random 
matrix model Q . We choose chiral representation 
of the Dirac gamma-matrices. Then the matrix 
of the Dirac operator in a certain basis has the 
block-diagonal form: 
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Nonzero entries X are approximated by a com- 
plex random NxN matrix with a Gaussian distri- 
bution: P(X) ~ exp{-iV Tr XX^}. To adapt 
the model to our case we add a chemical poten- 
tial: /U7o- 

To find G one calculates the quenched free en- 
ergy: 



V = -(lndetO -£>)), 
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and differentiates it with respect to the bare 
quark mass z: G — —dV/dz. Then p is given 
by: 
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Figure 1. The evolution of cuts of the function 
G{z) = (V'V') with increasing p from (Q). 

A standard way to deal with the logarithm in- 
side the averaging in (0) is the 'replica trick'. 
One does the calculation for arbitrary number of 
fermion species n: 

V n = - - ln(dct"(z - D)), (5) 
n 

and then takes n — > hoping that the limit 
is smooth. Doing this calculation we find that, 
first of all, G is a holomorphic function of z. 
It is given by a solution of a cubic equation: 
G(0 + G) 2 - p 2 ) - {z + G) = 0. It is analytic 
except for branch points, where 2 solutions co- 
incide, connected by cuts. The cuts follow lines 
where the value of Re [G 2 — ln((z + G) 2 — p?)} on 
two solutions coincide. The physical sheet is de- 
termined by an asymptotic condition at z — > oo: 
G — > 1/z, from the normalization of p. 

The singularities of G are shown in Fig. ^. 
For p? < 1/8 the cut is on the imaginary axis. 
At p? > 1/8 the branch points bifurcate and the 
cuts go into the complex plane. Above a certain 
value of pi — 0.278 . . . the cut no longer goes 
through z = — there is no discontinuity across 
the imaginary axis, i.e., the chiral symmetry is 
restored. 

4. n — » and quenched theory 

One can now compare this replica trick calcu- 
lation to the actual distribution of eigenvalues of 
the random matrix D in the quenched theory. It 
is obtained numerically, by plotting the eigenval- 
ues for an ensemble of random matrices (^) on 
Fig. ||. First of all, the eigenvalues do not stay 
on the imaginary axis for any nonzero pi They 
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spread in a blob of the x-width growing ~ /x 2 at 
small fi. Exactly the same behavior is observed 
in the quenched QCD calculations [0 and leads 
to the puzzling conclusion that /i c w and 
vanishes in the chiral limit m = because there 
is no discontinuity in G across the imaginary axis 
when /i > 0. 

Now, however, we are able to find the explicit 
solution of the model and see the nature of this 
quenched blob. First of all, we must learn that 
the limit n — ► of QCD with n quarks is not 
smooth and the quenched theory is not a good 
approximation to the full theory at ^ > 0. Sec- 
ond, we want to know what theory has quenched 
QCD as its n — ► limit. We find that it is a 
theory where each quark has a partner with the 
conjugate Dirac operator. The free energy reads: 
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ln(det"(z-L>)(z* - D*)). 
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This statement can be made rigorous even beyond 
the random matrix approximation using equa- 
tions (g) and: d 2 (z-X l ) = -(l/47r)Aln \z - X t \ 2 . 

Naively, when n — ► the conjugate quarks de- 
couple: V n , n -> V{z) + V*(z*) 1 and G = -dV/dz 
is holomorphic. In fact, the quarks ip and con- 
jugate quarks \ mrx f° r some values of param- 
eters z and n, i.e., precisely inside the blob. A 
condensate (xV 7 ) develops nonzero vacuum expec- 
tation value. This condensate carries nonzero 
baryon number. It was indeed observed in a 
Monte Carlo simulation of the SU (2) theory with 
quarks Q, which are self-conjugate in that case. 
We learn now that similar phenomenon occurs in 
the quenched theory. Calculating (o) we find § 
for the density of eigenvalues inside of the blob: 
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which one can check is in perfect agreement with 
the numerical data Fig. |[ 

5. Conclusions 

Using a simple random matrix model of chi- 
ral symmetry breaking we demonstrate a general 
phenomenon that at nonzero /x the limit n — > of 
QCD with n quarks is not smooth and does not 
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Figure 2. Eigenvalues of 20 random 100 x 100 
matrices on the complex plane z at two values of 
y 2 : (a) 0.06, (b) 0.40. The line is the boundary 
of the blob calculated analytically from @. 

coincide with the quenched theory. On the other 
hand, the quenched theory is a smooth limit of 
a theory where each quark has a conjugate part- 
ner. Another way of putting this is to notice that 
the difference between such a theory and QCD 
with In quarks is that the phase of the fermion 
determinant is omitted in the former. We learn 
that the phase of the determinant is extremely 
important for a simulation of full QCD at finite \x. 
Without this phase the result /x c = is natural, as 
happens in quenched QCD. This is due to the for- 
mation of the baryonic condensate, which breaks 
the (replica) chiral symmetry of quarks and con- 
jugate quarks. The Goldstone boson associated 
with this breaking is the so-called baryonic pion 
— a bound state of a quark and a conjugate anti- 
quark. The presence of this particle explains the 
result /i c = 771^/2. 
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